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1 Introduction 



The two-by-two ABCD matrix plays a central role in optical sciences. The 
four elements of this matrix are all real, and its determinant is one. Thus, it 
has three independent parameters. Yet, it has a rich mathematical content 
which could lead interesting results in physics. 

It is generally assumed that this matrix can be diagonalized by a rota- 
tion, but this is not the case as shown in our previous paper [TJ. We have 
shown there that this matrix can be brought to an equi-diagonal matrix by a 
rotation, and then by a squeeze to one of the following four Wigner matrices. 



cos 9 — sin6>\ /cosh A sinhA\ /l — -y\ /l 0\ 
sinfl cosfl J' UinhA cosh A/' \0 1 J' V7 l) ' ( j 



This squeeze portion of the similarity transformation is not yet widely known. 
Thus the similarity transformation which brings the ABCD matrix to one 
of the Wigner matrices is a rotation followed by a squeeze. Even though the 
two triangular matrices in Eq.(TfJ) can be similarity-transformed from each 
other, it is convenient to work with the four branches of the ABCD matrix. 

The purpose of this paper is to reduce these four matrices into one analytic 
matrix with four different branches. First of all, each of the matrices in Eq.flTJ 
is generated by 



1/0 -z\ 1/0 i\ 1/0 -i\ 1/0 

2 \i J ' 2 \i ) ' 2 V0 J ' 2 \i 



(2) 



respectively. The last two matrices can be obtained from a linear combina- 
tion of the first two, with two independent coefficients. We can then study 
the general property of the ABCD matrix by exponentiating the linear com- 
bination of the two matrices 

° "o ! ) - (! J). « 

and making Taylor expansions. 

One of the present authors noted this aspect of the ABCD matrix while 
studying optical activities [2] . He then concluded that the asymmetric optical 
activity can lead to the study of the fundamental space-time symmetries of 
elementary particles [31 H]. 

In this paper, we study the resulting exponential form more systemati- 
cally. We first exponentiate the linear combination of these two independent 
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matrices. While the exponent is a fully analytic function, the Taylor expan- 
sion of the exponential form leads to complications, leading to four separate 
branches. 

Again in this paper, we use the same optical activity to study the ori- 
gin of this branching property. We note then that the exponential form is 
convenient for repeated applications of the ABCD matrix, such as periodic 
systems including laser cavities and multi-layer optics. 

In Sec. [2j we discuss how the ABCD matrix can be written as an expo- 
nential function of one analytic matrix, with four branches. In Sec. [3j we use 
optical activities to study the physics of the mathematics of Sec. [2j Section 0] 
is devoted to application of this methods to periodic systems. Laser cavities 
and multilayer optics are discussed in detail. 

2 Exponential Form and Branches 

Let us start with the ABCD matrix as a rotated equi-diagonal abed matrix: 



[ABCD] = R(a)[abcd]R(-a) 



(4) 



where R(a) is a rotation matrix 





and [abed] is an equi-diagonal matrix 




(6) 



with 



tana = 



a = d 



D-A 
B + C 
A + B 
2 



b = 



C) + yf(A-D) 2 + (B + C) 2 
2 



c = 





2 
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Since the determinant of the ABCD matrix is assumed to be one, this 
matrix has three independent parameters. One of those parameters is the an- 
gle a. Thus, the abed matrix has two independent parameters. Since the two 
diagonal elements of the abed matrix are the same, it can be exponentiated 

as 

[abed] = exp{rM(0)}, (8) 

with 

— cos 8 + sin 8 ' 



M(8) = ( I — " ' . (9) 

V cos 8 + sm d / 

Here the two independent parameters are r and 8. Thus, we are led to study 
in detail this M{8) matrix which can also be written as 

M(0) = (° 1 " 1 )cos^+(J J) sin^ (10) 

Other than the factor of (z/2), this expression becomes the four generators 
given in Eq.(J2]) when 8 = 0, tt/2, tt/4, — 7r/4 respectively. 

In this way, we can combine the four Wigner matrices into an exponen- 
tial function of one analytic matrix. The problem is how to compute the 
exponential form of Eq.flH]). Its Taylor expansion is 



[^cd]=J2-r[M(8)} n - (11) 



This is an infinite series except at 8 = ±7r/4. If 8 = 45°, the M matrix 
becomes 

"-(AS)- w> 

Since M 2 = 0, the series truncates. The abed matrix becomes 

[abed} = (^ J). (13) 



Likewise, when 8 = — 7r/4, 



[oM = (J T f)- (14) 

This aspect is illustrated in Fig. [TJ The Taylor series truncates at 8 = 
±7r/4. In the circular regions 1 and 2, (sin^) 2 is smaller than (cos#) 2 . In the 
hyperbolic regions 1 and 2, (cos^) 2 is smaller than (sin#) 2 . 
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Figure 1: Forms of the ABCD matrix depending on the angle 9. 

Also in Fig. dj the symmetry of trigonometry tells us the region 7r/2 < 
9 < 37r/2 is redundant if we allow both positive and negative values of r in 
Eq.flH]). In this way, we restrict cos# to positive values. In circular region 1, 
sin 9 can be both positive or negative. 

In the region (cos9) 2 > (sin^) 2 , and 16*1 < 7r/4, if sin^ is positive, we can 
write the M matrix as 



M(0) = ^(20) ( ^ -«P(-»),. (15 ) 



with 



/cos^-sin^ 

exp [-T]) = J , 16 

V cos 9 + sin 9 
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where r] is positive. This formula is valid also when sin# is negative, but 77 
is also negative. We now write M{9) as 

M (»)=\^)( e T e °0(? VHT e-°/0- (17) 

Then 

r"/ 2 \ /0 -l\ n /e"/ 2 



{M{6)T = (cos(29))^ ( 6 "J e , /2 )(5 ^ e _° /2 ). (18) 



Thus 



which is 



, e"^ 2 \ /cos0 -sin^We"/ 2 , 
[°^ = ( e^Hsinc* coed A e^ 2 ) ' (19) 



[«h= ~ T ) (20) 

V e ' sin cos <p J 

with 

= r ^/cos(20). (21) 
In terms of the four parameters of the ABCD matrix, 

, ^ + 5 

COS0 = , 



e -2 V = ± = C-B-y/(B + C)* + (A-D)* 
c C-B + J(B + C) 2 + (A-Dy' 



where cos <fi is smaller than one, and b is negative. 

If (sin#) 2 > (cos#) 2 , or tt/A < \9\ < tt/2, we have to consider two separate 
regions in Fig. [U where cos 9 is positive, while sin 9 can take both positive 
and negative signs. cos(26 l ) is negative. The M matrix should becomes 

exp (—77) 



with 



M W = V-cos (20) ( ^ '»). (23) 



, . / sin - cos 6* 

exp(-r7) = 4 / — -, (24) 

V sm 9 + cos 9 

for both positive and negative values sin 9, but 77 is positive and is negative 
respectively. 
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Then the Taylor expansion leads to 
[abed] = 

with 



cosh x e v sinh \ 
e v sinh \ cosh % 



\ = / y-cos(20). 
In terms of the parameters of the original ABCD matrix, 

A + B 
coshx = — - — , 



(25) 
(26) 



-2 V 



b B-C + J(B + C) 2 + (A-Dy 



C-B + J(B + C) 2 + (A-Dy 



(27) 



Here cosh \ is greater than one, and both b and c are positive. 

We can now go back to Eq. ffTB"]) and Eq. (I24[) . and write tan# in terms of 
i]. Then tan# becomes 



tan# 



b + c (A-D) 2 + (B + C) 2 



c-b C-B 
for all values of 9 between —ir/2 and 7r/2. The parameter r is 

6 2 + c 2ll/2 



(28) 



-26c 

b 2 + c 2 
2bc 



1/2 



(29) 



for (sin#) 2 < (cos^) 2 and (sin^) 2 > (cos^) 2 respectively, with 



~b 2 + c 2 ' 


1/2 


~2{B 2 + C 2 ) + (A — D) 2 ' 


1/2 


2bc 




ABC + {A- D) 2 





(30) 



Let us now look at how the transition of the abed from Eq. (|20|) to Eq. (|25|) . 
This is a puzzling question because the matrix M{9) remains analytic in the 
neighborhood of 9 = it /A (see Fig. [TJ In order to tackle this problem, we 
write M{9) of Eq.© as 



M{9) 



COS I 





1 + tan 9 



1 — tan i 




(31) 
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In the neighborhood of 9 = 7r/4, we can set cos 9 = l/y/2 and (l + tan#) = 2, 
and 

M W = -L(° -(!-*»»)). (32) 
Then up to r 2 , the Taylor leads to 

If 9 is smaller than tt/4, the diagonal elements of this matrix are smaller 
than 1, like cos<ft in Eq. (1201) . If 9 becomes greater than 7r/4, the diagonal 
element becomes greater than 1 like coshx in Eq.( l25l) . If tan# = 1, the result 
becomes that of Eq. ( TTBl . 

We can give a similar reasoning for the neighborhood of tan# = — 1. The 
Taylor expansion leads to 

r h A\ - n-r 2 (l + tan^)/2 -rV2 \ 

[aDCai ~ \ r(l + tan^)/ v / 2 1 - r 2 (l + tan0)/2 ) ' [ > 

leading to Eq.dHJ) for 9 = -7r/4. 

The exponential form given in Eq.® is very convenient when we study 
periodic systems where the ABCD matrix is applied repeatedly. We shall 
return to this problem in Sec. HI 



3 Optical Activities 

In his recent paper [2] , one of the present authors used the two- by-two matrix 
formulation of optical activities applicable to the transverse electric field of an 
optical wave. The direction of the electric component rotates as the optical 
wave propagates. In the real world, the medium causes also an attenuation 
of the transverse components. This does not interfere with the rotational 
character. However, there is a problem if the dissipation coefficients are 
different for two perpendicular directions. 

Let us start from a circularly polarized light wave which can be decom- 
posed into the right-polarized and left polarized components. If they have 
different indexes of refraction, we can write the light wave as 

(#0 = 2 (!) GXP ^ ( klZ ~ + 2 (-i) ex P0(^-^*)} (35) 
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This two terms can be combined into 



with 

k= 1 -(k l + k 2 ), 7 = ^^. (37) 
If we start with a polarized light wave taking the form 



E x \ _ ( Aexp {i(kz — oot)} 



EyJ V 



(38) 



the optical activity is carried out by the rotation matrix 



\ sin(7z) cos(72) J ' 

The optical ray is expected to be attenuated due to absorption by the 
medium. The attenuation coefficient in one transverse direction could be 
different from the coefficient along the other direction. Thus, if the rate of 
attenuation along the x direction is different from that along y axis, this 
asymmetric attenuation can be described by 

exp {-Hiz) \ _ Xz f exp faz) \ , , 

exp(-fi 2 z)J \ exp(-fiz))> 

with 

/i 2 + /ii /i 2 -/ii (A . 

A = ^— , »=—2— ( 41 ) 

The exponential factor exp (—Xz) is for the overall attenuation, and the ma- 
trix 

' exp (j^z) 
exp (—fj,z) 

performs a squeeze transformation. This matrix expands the x component 
of the polarization, while contracting the y component. We shall call this 
the squeeze along the x direction. 

The squeeze does not have to be along the x and ^/directions For con- 
venience, let us rotate the squeeze axis by 45°. Then the squeeze matrix 
becomes 

S^z)=( COS ^ Z ] Sin ;HV (43) 
V smh(fiz) cosh(fiz) ) 



(42) 
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If this squeeze is followed by the rotation of Eq. (l59"|) . the net effect is 



Xz ( cos(72) — sin(72;) \ / cosh(/z2:) smh(/iz) \ 



sin(7z) cos(72;) / \ sinh(//z) cosh(/x2;) J 



(44) 



where z is in a macroscopic scale, perhaps measured in centimeters. However, 
this is not an accurate description of the optical process. 

This happens in a microscopic scale of z/N, and becomes accumulated 
into the macroscopic scale of z after the N repetitions, where N is a very 
large number. We are thus led to the transformation matrix of the form 



-Xz 



/N S{^z/N)R{ 1Z /N) 



N 



In the limit of large N, this quantity becomes 



-Xz 



( 1 »z/N\ 
\fjiz/N 1 / \^z/N 



-7Z/JVN 



N 



Since ^z/N and fiz/N are very small 
Z(j,fi, z) -- 

For large N, we can write this matrix as [2] 



-Xz 



1 
1 





1 + v 



-7 + M _?L 

N 



i N 



Z(j, n, z) = e~ Xz exp {kzM(6)} 



where the M matrix is 
M{9) 

with 





cos 9 + sin 9 



cos 9 + sin 9 




k = 
cos 9 

sin 9 



7 2 + /i 2 , 



7 



(45) 



(46) 



(47) 



(48) 



(49) 



(50) 
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We note here that the M(6) matrix of Eg. (149 p is the same as that of Eq.flH]) 
which determines the branch property of the ABCD matrix. At this point, 
it is more convenient to work with kM(6). 

If 7 > fi, the tM matrix can then be written as 

kM= v¥ r 7(° ~l~ V ), (52) 
where rj of Eq. ffl6|) becomes 



e - = = . C03 ^- Si ^ . (53) 

Y 7 + /i V cos 6 1 + sin 6> 



Thus, the exponential function in Eq. (l48j) can be evaluated according to the 
procedure defined in Sec. EJ This expression is the same as that of Eq. (116j) . 
The exponential form exp (kzM) in of Eq. (|48p becomes 



/ cos^'z) — sxn^'z) \ ( e v ^ 2 
e^ 2 J V sin(y^) cos(Vz) ) \ e~^ 2 J ' 



(54) 



where 



V = yV - /i 2 - (55) 
The transformation matrix of Eq. fj48|) takes the form 



l7 '^' j U"sin(7'z) cos(Vz) ' ' 1 ' 

If fi > 7, the rM matrix becomes 



where 77 of Eq. (124]) takes the form 



-2 V _ _ / sing-cosfl 

e - VJTR" Vcos^ + sin^' (5bj 
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and Z becomes 



Zh a z) - e- Xz ( cosh (^) e "smh(^)\ , . 

l7 '^' > U" sinh(/x'z) cosh(jz'z) /' 1 J 

where 

= \/V 2 -7 2 - ( 6 °) 

In this section, we discussed a system of optical activities with asymmetric 
dissipation as a physical illustration of the mathematical procedure discussed 
in Sec. [2] We have already seen that M(9) of Eg. (149D has the same form as 
that of Eq.flH]), and that angle 9 can be defined in terms of the parameters 7 
and n, as shown in Eq.f l50p . The parameter rj can also be defined in terms 
of 7 and (jl, and its expression is the same as the one given in terms of the 
angle 9. 

As for the branches, we note that both 7 and \x can be negative and 
positive. Thus, the angle 9 can cover the entire range from zero to 2tt. We 
can write 7' and pi as 



7' = k\J cos 2 9 - sin 2 9, // = k\J sin 2 9 - cos 2 9. (61) 

Since cos 2 9 — sin 2 9 = cos(26 l ), j'z and fi'z correspond to <fi and x °f Eq. fl21l) 
and Eq. (l26l) respectively. 

If we start with \x = 0, it is simply a rotation of the transverse component 
of the electric field and the overall attenuation factor is exp (— Xz). The rate 
of this rotation decreased as \x increases, and the rotation stops at 7 = fi. 
For fi > 7, there are no rotations. It would be very interesting to test these 
effects experimentally. 

We should not forget the fact that the equi-diagonal [abed] matrix is a 
rotated ABCD matrix. The rotation matrix is given in Eq. (jSJ). This rotation 
changes the optical ray of Eq.(l36i) to 

This is also an observable effect. 

We have seen in this section that the asymmetric optical activity can 
serve as an analog computer for the mathematical procedure given in Sec. [2] 
which is in fact an alternative to the diagonalization of the ABCD matrix. 
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4 Periodic Systems in Optics 

Let us summarize what we can do about the ABCD matrix. 

1. We should first rootate to an equi-diagonal matrix [abed]. 

2. If the diagonal elements of this equi-diagonal matrix are smaller than 
one, it can be written as 



with exp (—rj) = —b/c, which can also be written in terms of the ele- 
ments of the original ABCD matrix, as shown in Eq. fl2~Tj) . 

3. If the diagonal elements of the equi-diagonal matrix are greater than 
than one, the matrix can be written as 



with exp (—77) = b/c, which takes the form of Eq.f l26l) in terms of the 
elements of the ABCD Matrix. 

4. If one of the off-diagonal elements vanish, the diagonal elements have 
to be one. 

5. It is possible to combine all these cases into one exponential function 
of one analytic matrix. It can be written as 




cos (j) —e v sin 
v sin <p cos <p 



(63) 




(64) 



[abed] = exp < r 



( 





cos 9 + sin 9 



— cos 6 + sin 9 




(65) 
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When 9 = ±7r/4, the Taylor series truncates, and 



[abed] 



( 



1 



r(l T l)/ v / 2 



) 



(66) 



r(l±l)A/2 



1 
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Figure (a) 



The cycle can start anywhere 
between the lenses. 

r 




The cycle starts at the midway between 
the lenses when [ABCD] is equi-diagonal 



Figure (b) 



a> 
o 



midway 



.(S/2). 



o 

o 



(s-x) 



Figure 2: Optical rays in a laser cavity, (a) Multiple cycles in a laser cavity 
are equivalent to the beam going through multiple lenses, for which one cavity 
cycle corresponds to the propagation of light through a sub-system of two 
lenses. The ABCD matrix becomes equi-diagonal when the cycle begins at 
the midway between the lenses, (b) A laser cavity consisting of two concave 
mirrors with separation s. 



4.1 Laser Cavities 

A laser cavity consists of two concave mirrors separated by distance s as 
illustrated in Fig. [2j The mirror matrix takes the form 



1 

-2/R 1 



(67) 



where R is the radius of the concave mirror. The separation matrix is 

J 0- («> 
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If we start the cycle from one of the two mirrors one complete cycle consists 

If we start the beam at the position x from the mirror, then one complete 
cycle becomes 

1 iW 1 0\(ls-x 
I J \-2/R l)\0 1 

1 x\ ( 1 Wl s-x\ 
l)(-2/JJ l)U 1 )■ (70) 

This cycle consist of two identical half cycles. 

Thus, we shall use the half-cycle matrix as our starting point. Then the 
half-cycle ABCD matrix becomes 

i^i = (J i)(- 2 % ?)G '!*)■ < 71 > 

It is now possible to replace replace R and x and by Rj s and x/s respectively 
and set s = 1 pp. Then 

[ABCn^^-Jf 1 -_%'$:$). (72) 

where / = s/R, and is expected to be a small number because the mirror 
radius R is much larger than the separation of the mirrors. 

It can be brought to an equi-diagonal form by a rotation as given in 
Eq.fllj). According to Eq.([7j), the rotation angle is 

2f(2x - 1) 

tana = ^ J —^. 73 

l-2f(l + x-x 2 ) K J 

This angle is zero when x = 1/2. In this case, the laser cycle starts at the 
midway between the lenses pp. Then the ABCD matrix becomes 

[aM = ( 1 _a/ Y- 7 / 2 )' (74) 



This matrix can then be written as 



■ —e 1 sin d) cos <t> ' 
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with 



cos0 = l-/, e2V = 2 -Jf- ( 76 ) 

This is the result we obtained in our earlier paper on laser cavities [5] , where 
the cycle starts from the midway between the lenses. The signs of <fi and 
r) are opposite to those given in Eq. fl20|) . but this is purely for convenience. 
There are no fundamental problems. 

We can now write this expression in an exponential form 

W =ap f(-J +B hi« cos 61 + am 61 \\ (7?) 



with 



tai ^ = 2T3f 
13/ - 17/ 2 



8 - 17/ 2 



1/2 

0, (78) 



where <p is given in Eq. f[76|) . Since the radius of the mirror is much larger 
than the mirror separation, / is a small number, and tan# is close to one 
and r is a small number. 

The TV-cycle laser consists of 2N half-cycles, and its abed matrix is 

[ M- = exp{ 2JVr (_ sin ; +cos() -«;*')}. (79) 

This section is a straight-forward application of the procedure given in 
Sec. [2j We know that sm(r8) 2N is not sin(2Nr6) , but the exponential form 
gives us the convenience of [exp (ir8)] 2N = exp (J,2Nr9). We have given a 
two-by-two matrix formulation of this convenience applicable to the ABCD 
matrix. 



4.2 Multilayer Optics 

From the physical concept of Wigner's little group whose transformations 
leave the four-momentum of a given particle invariant [31 H], it has been 
established in the literature that [6] 

S( V )WS(- V ) = R(0B(-2X)R(0, (80) 
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where W is one of the four Wigner matrices given in Eq.flT]), R(£) is the 
rotation matrix of the form of Eq.([5]), and 



B(X) 









-n/2 



cosh(A/2) sinh(A/2) \ 
sinh(A/2) cosh(A/2) J 



(81) 



and the continuous parameters £ and A take care of the four different Wigner 
parameters. These parameters can be written in terms of r\ and the parameter 
of the Wigner matrix, as shown in Ref . [TJ [6] . 

Since the left side of Eq.f l80p can be written as an exponential form, we 
can write 





, cos 9 + sin 9 
where r and 9 are also continuous variables. 

I 

Beam 



— cos 9 + sin 9 




(82) 



Medim 1 



J 



Cycle starts at the boundary 
between medium 2 and medium 1. 



Figure 3: Multilayer consisting of two different refractive indices. One com- 
plete cycle starts at the boundary between medium 2 and medium 1. 

With this mathematical preparation, let us study multilayer optics. In 
this branch of optics, we have to consider the ABCD matrix applicable to 
two beams moving in opposite directions, one which is the incident beam 
and the other is the reflected beam [7]. We can represent them as a two 
component column matrix 



E_ 



—ikz 



(83) 



where the upper and lower components correspond to the incoming and re- 
flected beams respectively. For a given frequency, the wave number depends 
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on the index of the refraction. Thus, if the beam travels along the distance 
d, the column matrix should be multiplied by the two-by-two matrix [7] 

P(Pl) = { o e -l/2), (84) 

where /3j/2 = kjd and j is denoting each different medium. If the beam prop- 
agates along the first medium and meets the boundary at the second medium, 
it will be partially reflected and partially transmitted. The boundary matrix 
is0 

B(u) = ( C0S ^ Si t ( r^)' ( 85 ) 

Vsmh(z//2) cosh(z//2) / 

with 

cosh(z//2) = l/t 12 , sinh(z//2) = r 12 /t 12 , (86) 

where £i 2 and r 12 are the transmission and reflection coefficients respectively, 
and they satisfy (rf 2 + £^ 2 ) = 1. The boundary matrix for the second to first 
medium is the inverse of the above matrix. Therefore, one complete cycle, 
starting from the second medium, consists of 

B(v)P{p x )B(-v)P(fa), (87) 

as illustrated in Fig. [3j This complex-valued matrix can be cast into a real 
matrix by a similarity transformation with the transformation matrix 

1 / e i7r / 4 e^/ 4 

This transforms the boundary matrix B{y) of Eq.( )85|) to a squeeze matrix 
S(v) of Eq. llHT]) . and the phase shift matrices P((3j) of Eq. flM]) to rotation 
matrices R(f3j) of the form given in Eq.([5]). We are thus led to consider the 
ABCD matrix of the form 

[ABCD] = S{u)R{f3 1 )S{-u)R{f3 2 ). (89) 

If W in Eq.f lSU]) is a rotation matrix, we can write 

S{v)R{^)S{-v) = R^ X )B{-2\)R{^) (90) 

where 



cosh A = (cosh i/) \jl - cos 2 (/3i/2)tanh 2 ^, (91) 

c cos(/3i) 

cos 6 = 1 ■ (92) 

(coshz/Wl — cos 2 (/9i/2) tanh v 
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The ABCD matrix can then be simplified to 

[ABCD] = Rfa)B{-2\)Rfa) (93) 

with 

£2 = 6 + h (94) 
It is now possible to write the above form as 

[ABCD] = R(a)[R(£)B(-2\)R(£)]R(-a), (95) 

with 

£ = ^(6 + 6), « = ^(£i-6). 

The role of the rotation matrix R(a) matrix is clearly defined in Sec. El Thus 
R(£)B(—2\)R(£) is the equi-diagonal matrix, and 



R(0B(-2\)R(0 



cosh A cos £ — (sin £ cosh A + sinh A) 

sin £ cosh A — sinh A cosh A cos £ 

(96) 

Thus, if (cosh A cos £) is smaller than one, we can write this matrix as 



cos <p — e 1 sin 
~ v sin 6 cos 6 



■nu 

(97) 



with 



2v (cosh A) sin £ + sinh A 



cos0 = (cosh A) cos £, e v = — — — -. (98) 

(cosh A) sin 4 — smh A 

Thus, if (cosh A cos £) is greater than one, we should write the equi-diagonal 
matrix as 

cosh x —e v sinh \ ' 
-e~ v sinh \ cosh \ 



(99) 



with 



2v sinh A + (cosh A) sin £ 



cosh x = (cosh A) cos f, e v = — — — — (100) 

smh A — (cosh A) sin £ 

We are now interested in the exponential form 
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with 



As for the r parameter, 



tan 6* 



tanh A 



sin 2 £ + tanh 2 A 
sin 2 £ — tanh 2 A 



1/2 



sin 2 £ + tanh 2 A 
tanh 2 A — sin 2 £ 



1/2 



X, 



(102) 



(103) 



for (cosh A cos £) < 1, and for (cosh A cos £) > 1, respectively. 

In this section, we started with two media with two different indexes of 
refraction, corresponding to two rotation matrices R((3i) and R{^2) given 
in Eq. (l89j) respectively. However, the combined effect in not necessarily a 
rotation matrix. It can be analytically continued to the hyperbolic branch 
through the exponent of the ABCD matrix. 

When (cosh A cos £) 2 = 1, one of the off-diagonal elements in Eq.f llOip 
vanishes, and this case was repeatedly discussed in the literature [TJ [6], also 
in the present paper. 



Concluding Remarks 

In this paper, we noted first that the two-by-two ABCD matrix can be 
represented as a similarity transformation of one of the four matrices which 
we choose to call the Wigner matrices. We then combined these Wigner 
matrices into one exponential form of an analytic matrix. 

While cos0 and coshx correspond to a circle and a hyperbola respec- 
tively, the lines in Fig. [1] correspond parabolas in the four-dimensional repre- 
sentation of the Lorentz group [21 E]- Ancient Greeks used a circular cone to 
combine these curves into one. This is the reason why we call them conic sec- 
tions. It is gratifying to note that the optical devices we discussed in Sees. [3] 
and H] can play the role of a conic section. Instead of three-dimensional cone, 
we used a two-dimensional plane in Fig. [TJ 

We have seen in this paper that Taylor expansion of this analytic form 
results in four branches. How does this happen? Let us go to the Taylor 
expansion of Eq. ffTT]) . This infinite series truncates at (sin#) 2 = (cos^) 2 
or along the two lines in Fig. [TJ We are not familiar with mathematical 
singularities resulting from the truncation of the infinite Taylor series. This 
appears to be an interesting problem in mathematics, but it is beyond the 
scope of this paper. 
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